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We consider classes of translationally invariant black hole solutions whose equations of state 
closely resemble that of QCD at zero chemical potential. We use these backgrounds to compute 
the ratio (/s of bulk viscosity to entropy density. For a class of black holes that exhibits a first 
order transition, we observe a sharp rise in ^/s near Tc- For constructions that exhibit a smooth 
cross-over, like QCD does, the rise in C^/s is more modest. We conjecture that divergences in C^/s 
for black hole horizons are related to extrema of the entropy density as a function of temperature. 
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The anti-de Sitter / conformal field theory (AdS/CFT) 
correspondence [E 13, 9 has generated interest in using 
thermal A/" = 4 super- Yang-Mills theory (SYM) to under- 
stand quantum chromodynamics (QCD) at finite tem- 
perature. A conspicuous shortcoming of this approach 
is precisely the conformal invariance of SYM. This im- 
plies, for instance, that in the SYM plasma, the speed 
of sound Cs equals l/-\/3 and that the bulk viscosity C 
vanishes at all temperatures. QCD only exhibits confor- 
mal behavior in the high-temperature regime. In order 
to roughly capture the behavior of QCD across a larger 
range of temperatures, we are led to consider gravity du- 
als of gauge theories that break conformal invariance [3l| . 
The minimal action on the gravity side that can describe 
Lorentz- invariant, non-conformal theories is 



S=^J d'x^(^R-^{dq>r-Vi^)y (1) 

where V{(j)) is the potential for the bulk scalar field (j), 
and = y/SnG^ is the five-dimcnsional gravitational 
constant. We restrict our attention to backgrounds of 
the type 



ds^ = e^^^^^ {~h{r)dt^ + dx"") + e2^('') 
(/) = (/)(r) dx^ = {dx^f + {dx^f ^ 



dr^ 

W) 

- {dx^f . 



(2) 
(3) 



This form is the most general ansatz with translational 
symmetry in the {t, x) coordinates and 50(3) symmetry 
in the x directions, as is appropriate to describe an infi- 
nite static thermal plasma. The equations of motion for 
the functions A, B, and h come from plugging the ansatz 
(l2|)-([3]) into the equations of motion following from the 
action ([1]). See in which techniques for finding solu- 
tions of the form ©-([S]) and the corresponding equation 
of state are more fully explained. 



'Electronic address 
t Electronic address: 
Electronic address: 
§ Electronic address: 



ssgubser@Prince ton.EDU| 
anelloreOPrinceton. EDUT 
spufuOPrinceton.EDUl 
frochaQPrinceton.EDU] 



The AdS's-Schwarzschild solution can be recovered as 
the limit of the above construction where cj) vanishes and 
V{4)) is a negative cosmological constant term. More 
generally, if at small (j) one has 



n0) = -^ + 
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then the gravity solution ([I])-® will be asymptotic to 
anti-de Sitter space with radius L. An asymptotically 
AdS spacetime on the gravity side is equivalent to con- 
formal invariance of the field theory in the UV. Gravity 
backgrounds constructed from potentials which satisfy 
(|4|) are dual to relevant deformations of a conformal field 
theory: 



£ = £ 



CFT 



(5) 



where A^ is the energy scale of the deformation and A is 
the dimension of the operator dual to (f>. According 
to the AdS/CFT dictionary, A can be identified with the 
larger root of 



A(A - 4) 



2r2 



(6) 



We will only be interested in the case 2 < A < 4, 
which corresponds to relevant deformations that obey the 
Breitenlohner-Freedman (BF) bound 0, 0, 0|- 

A background of the form (H])"© has an event horizon 
if h has a zero. Let th be the value of r closest to the 
conformal boundary where h vanishes. Thermodynamic 
quantities such as entropy density s and temperature T 
are parameterized by rn'- 



T : 



47r 



The speed of sound Cg can be computed from 
" dlogs 



(7) 



(8) 



We exclude from consideration nonzero chemical poten- 
tial for baryon number. To include this, we would have to 
add a gauge field to the action ([T]) and consider charged 
black holes. 



2 



liV = Vqc'''^ with Vq < 0, then the equations of motion 
following from ([1]) can be solved analytically ^ , and the 

speed of sound is constant: = ^ — But these black 
holes aren't asymptotically anti-de Sitter because V{(j)) 
has no maximum. If instead V{4') interpolates smoothly 
between ^ for small (j) and Voe"''^ for large cj), then the 
black hole solutions have a temperature-dependent speed 
of sound, Cs{T). In the mapping between V{(j)) and 
Cs{T) is explored in some detail. Within certain limits, 
given Cs{T), one can find a V{(f)) to reproduce it using 
black holes. 

It probably isn't possible to obtain an arbitrary V{(j)) 
from string theory. However, it is typical in gauged su- 
pergravity to find potentials with local extrema and ex- 
ponential increase or decrease as canonically normalized 
scalars become large. In any case, it is our goal to design 
a potential V{(j)) to reproduce the equation of state of 
QCD. 

It is perhaps surprising that the simple potential 



-12 cosh 70 - 



(9) 



lagrangian 



C 



QCD 



fcrmionic terms . 



(10) 



where ao = g^/^ir is the bare strong coupling constant. 
To compute the dimension of tr F"^^ , one starts by notic- 
ing that the trace of the QCD stress-energy tensor 



^trF^ 



fcrmionic terms 



(11) 



is RG-invariant, so it should scale classically. In (fTTI) . a 
is the renormalized coupling at scale Q, and /3(a) is the 
QCD beta-function 



(12) 



For any operator O, dO / d\ogQ — —OA, where A is 
the sum of the classical and anomalous dimensions of 
O. Thus, differentiating pTjl with respect to logQ, one 
obtains 13211 



4 + /?'(«) 



2/3(a) 



(13) 



with 7 w 0.606 and h w 2.057 approximately reproduces 
the squared speed of sound versus temperature as derived 
from lattice data on 2 -I- 1-flavor QCD: see figure [TJ Be- 
cause the equation of state exhibits a cross-over rather 
than a sharp phase transition, we have to use some pre- 
scription to determine Tc in order to plot versus T /Tc- 
We define Tc as the inflection point of s/T^ as a function 
of T. 

The quoted value for 7 corresponds to setting = 0.15 
in the extreme IR. Although hadron resonancegas mod- 
els give values of ranging as high as 0.2 (see [o'l and ref- 
erences therein), the value cj. = 0.15 is in a phenomeno- 
logically interesting range. The equation of state follow- 
ing from ([9]) is fairly close to the quasiparticle model of 
, based on a chiral extrapolation of lattice data. The 
black hole model is complementary to a quasiparticle de- 
scription in that it should work well precisely when no 
weakly coupled quasiparticle description is available, re- 
minding us of the correspondence principle of [1G| . Thus, 
the picture we advocate in using the potential ^ is that 
the approximate validity of a black hole description of 
QCD is not lost suddenly during the smooth cross-over, 
but instead gradually, so that the black hole continues to 
give an approximate guide to the dynamics at least down 
to Tc, and perhaps even somewhat below it. This is a de- 
parture from a more traditional picture, inspired in part 
by large N counting, where there is a sharp transition 
(usually first order) between a black hole description of 
a deconfined phase and a horizon-free description of the 
confined phase: see for example [Tl| . 

The quoted value for h corresponds to setting A « 
3.93, which is the dimension of F^^ in QCD computed 
at three loops and energy scale Q = 3 GeV. Here, 
is the rescaled field strength that appears in the QCD 



Reference 12 1 contains the exact expressions for P{a) 
and a{Q) At Q ^ 3 GeV, we get a « 0.253 and 

A « 3.93 at three loops. 

In summary, the UV matching to QCD doesn't at- 
tempt to capture asymptotic freedom, which probably 
requires going beyond the supergravity approximation; 
instead we match onto QCD at a finite scale, above which 
asymptotic freedom is replaced by conformal invariance. 

Black holes constructed with the potential ([9]) may 
slightly underestimate the rapidity of the cross-over of 
QCD. We therefore consider an alternative potential, 

-12 COSh7</. + 520' + &40^ + ^60' ,,,, 

VW = j2 ' (14) 

with 7 « 0.606, 62 ~ 1.975, 64 ~ -0.030, and ba « 
—0.0004, where the specific values are chosen so as to 
sharpen the cross-over almost to a second order phase 
transition. One can see from figure [T] that cj for the 
choice (IT4t is a close match to pure glue data of [l3| for 
T > Tc, so its behavior close to Tc is probably sharper 
than QCD's. But, by design, it still has a behavior rem- 
iniscent of hadron gas phenomenology for T < Tc. 

The shear viscosity of all black hole solutions we con- 
struct satisfies rj/s — 1/Att [l^ because we exclude higher 
derivative terms from the action. This low value of -q re- 
minds us that the regime of validity of a black hole de- 
scription cannot extend too far above Tc or too far below 
it. The bulk viscosity can also be studied (see for exam- 
ple [H, [H, [13, [i&])j and it is particularly interestin g to 
inquire how it behaves near Tc. There is a proposal [l9| 
that QCD exhibits a sharp rise in C/s close to the decon- 
finement transition, signaling "soft statistical hadroniza- 
tion" of the QGP. See also the earlier works [l^, l2l| . 
which deal with pure glue. 
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Bulk viscosity can be computed from the Kubo formula 

(15) 

GflH= I d''xdte'^'e{t){[Tu{t,x),Tkk{0,0)]) , (16) 



^ ^ lim -lmGfl(w) 

9 Lj^o UJ 



where Gr^uj) is the retarded two-point function of the 
spatial trace Ta of the stress-energy tensor. Two-point 
functions of the stress-energy tensor can be computed 
within AdS/CFT by examining the metric perturbations 
Sgjj around the background using the recipe of 

[2 2( 1 and subsequently justified in ]23| starting from the 
more fundamental prescription of z, 3]. In our case, the 
relevant metric perturbations exhibit rotational symme- 
try and thus only consist of dgoo, Sga, (Jgss, Sg^^, and 54>. 
Without loss of generality, we assume 5gii = 5g22 — Sg^^. 
Henceforth we work in the gauge r = cj), which is conve- 
nient because the equation for Sgu decouples from the 
other perturbations. Setting hn = e~^^Sgii, we find 
that the equation for Sgu reduces to 



1 



AA' + 3B' 



--2A+2B 



/l2 



h' 

6hA' 



h'B' 



hii, (17) 



where primes denote derivatives with respect to <j). We 
impose the normalization condition hn « 1 at the 
boundary of AdS, as well as infalling boundary condi- 
tions at the black hole horizon: 



hii 



\4> - 4>H 



(18) 



At any value of (j) between the horizon at = 4>h and 
the conformal boundary at = 0, the number flux of 
/ill quanta with frequency uj falling into the black hole 
is given by the conserved quantity 



^AA-Bi 



|Im/in/iii| 



(19) 



Given the number flux one can compute the imag- 

inary part of the retarded two-point function of Ta from 



Im Gb{^) 



(20) 



Heuristically, (pO|l is the statement that dissipation in the 
boundary theory is related to the probability for particles 
from the conformal boundary of AdS to be absorbed into 
the horizon. 

From (fT5|) and ([T9)l , one straightforwardly finds 



which gives 



-11 



4A'(0ff)2 



V 



OH) 



47r 



--11 1 



(21) 



(22) 



In deriving p2)) we have used the relation A'[<j)H) = 
-~V{4)H)/'iV'{4)H), which follows from the equations of 
motion for the background ©-([S]). Bulk viscosity mea- 
sures the hysteresis in nearly adiabatic S'0(3)-invariant 
perturbations of the thermal medium. Therefore, to ex- 
tract Cj we can compute the quantity c^^ appearing in 
([221) in the w ^ limit. 

A more detailed explanation of the derivation of (|22p . 
as well as a description of the numerical computation of 
Cii, is given in [23. Figure [1] includes a plot of C/s for 
the potentia ls pt and (fH)) . as well as for pure glue as 
obtained in [2l| from lattice simulations, and for QCD 
as obtained from the sum rule approach of [H, [13] . 

An alternative approach to constructing black holes 
that describe the thermal phase of non-conformal gauge 
theories has recently been suggested in [2^] , following ear- 
lier work [23, . The setup is similar to ([T]) , except that 



at large (j). It is argued that fluctu- 



ations around the zero-temperature background (which 
is singular) yield a glueball spectrum with ^ n for 
large excitation numbers n, suggestive of linear confine- 
ment. The black hole solutions of ^26j have a minimum 
temperature Tmin, where the specific heat diverges and 
vanishes. The equation of state resembles that of pure 
Yang-Mills theory, with a sharp transition at some tem- 
perature Tc- We will assume that Tc = Tmin, although it 
is more generic to have a first order transition at some 
Tc > Tmin- The potentials employed in [H [l^, [H] do 
not always have maxima, but the behavior of black holes 
close to Tr is similar to what one finds with 



y(0) = 



-12(l + a02)i/4(,Qgj^ U + 



L2 



(23) 



Potentials of this form also exhibit a glueball spectrum 
with ~ n. 

Because the phase transition is sharp for potentials of 
the form one might expect the behavior of C/s near 
Tc also to be sharp. And so it proves, as shown in figure[T] 
But C/s remains finite at Tc, even as the specific heat di- 
verges, similar to behavior reported in |18| . Exploration 
of potentials similar to (^5)) reveals the following related 
behaviors: 



• When the potential (|23p is modified to match the 
equation of state of pure glue more closely, the peak 
in C/s becomes broader and lower, and C/s becomes 
bigger well away from Tc- 

• It appears that (/s never diverges as long as > 0, 
which is equivalent to positive specific heat in the 
absence of chemical potentials. Potentials of the 
form exhibit solutions with < 0. In some 
cases, — > —00, corresponding to a minimum of 
the entropy density, and then (/s does diverge. 
Similar behavior arises when V{(j)) has a narrow 
region of sharp decrease. Of course, when < 0, 
the significance of ^ is more formal, because one is 
perturbing around an unstable background. 
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FIG. 1: LEFT: A comparison of different c^(r) curves. The solid red curve corresponds to the potential ([9| (Type I black holes). 
The dashed magenta curve corresponds to the potential (I14p (Type II black holes). The dot-dashed orange curve corresponds 
to the potential in (|23|) (Type III black holes) with a = 1 and b adjusted so that the dimension of the operator dual to <j!> is 
A « 3.93. We also show lattice results for pure glue (solid black curve) and 2 + 1-flavor QCD (open black squares), as well as 
Cs(r) for a 2 + 1-flavor quasiparticle model (QPM). The pure glue curve is based on [1^ and private communications from 
F. Karsch. The 2 + 1-flavor lattice QCD points are based on [2J|; for these points, we take Tc — 187 MeV, as estimated from 
the halfway point of the initial rise in an (e — 3p)/T* curve from [2^ ]. This differs from the value Tc — 196(3) MeV computed 
by a different method in The QPM points are based on 

RIGHT: A comparison of ("/s results for the black holes described above (see legend) and the results of [20l.[2H. Lattice results 
for pure glue from are shown in blue. The solid curves correspond to the sum rule result for QCD with 2-1-1 flavors from 
[20l |. using three representative values of the frequency parameter lvq. 



When V'{(j))/V{4>) is slowly varying, (/r] sa 2(i — 
Cg). This adiabatic approximation can be derived 
by dimensional compactification of a conformal 
field theory, and it is generally a good indicator 
of the order of magnitude of C, even when < 0. 
See [^,[i2[ for related observations, and in partic- 
ular [1^ for the conjecture that C/v ^ 2(i — c^) for 
all black hole solutions. 



In conclusion: five-dimensional gravity coupled to a 
single scalar contains the minimum amount of freedom 
needed to match an equation of state to a family of black 
holes. Constructions similar to those of [1^ show even 
sharper behavior in the equation of state near Tc than 
pure glue does. These constructions lead to a sharp rise 
in ^/s near Tc, though (^/s remains finite at Tc. The 
steepness of this rise is associated with the proximity 
of a minimum of s on a branch of thermodynamically 
unstable solutions. More generally, we conjecture that 
C/s diverges precisely when and where s has an extrcmum 
[S^ l . Constructions of Q , in which it is assumed that the 
smooth but rapid cross-over of QCD is described in terms 
of smooth but rapid cross-over behavior of dual black hole 
solutions, lead to a more modest rise in (/s near Tc. 

The steep rise oi (/s observed for the potential (1^51) 
is at least broadly consistent with the results of 0, l2ll |. 
although it is perhaps troubling that blunter behavior 
arises for potentials like that more closely match 



the pure glue equation of state. 

The more modest rise of C/s observed for the potential 
(ini) and its variant (|14p indicates some tension with the 
results of [l^. Our results suggest that (/s < 0.1 at Tc 
for QCD. If C/s is significantly bigger there, it either 
means that there is a subtlety in the matter lagrangian 
that we have not understood, or that Einstein gravity 
does not adequately describe the approach to confine- 
ment from above. If, instead, our results are closer to 
the true behavior of real-world QCD, it suggests t hat 



20| 
21| 



the parametrization of the spectral function in 
is somehow misleading, and that the lattice study 
needs to be extended to include fermions before being 
directly compared with QCD close to Tc. 
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